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1. INTRODUCTION 

The current decreases to a low value by increasing voltage to the maximum to retain the power 
constant. Line losses are also decreased in this manner. This means that less power is necessary for 
transmitting electricity via transformers at higher voltages. For this reason, power transformers are the most 
suitable choice. There are many types of transformers and transformers have many uses. But the basic 
purpose of transformers is to raise the voltage at the end of the generation and decrease the voltage at the end 
of the consumer [1], [2]. 

The mismatch between the transformer and its model is important specially during fault study. 
Faults are known as internal and through faults in the transformer. Internal faults are faults that occur within 
the transformer that can seriously damage the transformer's insulation and cause the transformer to break 
down. So, the transformer should be protected from these faults immediately. Such faults are split into 
electrical and mechanical faults [1]-[11]. The model for the transformer is of critical importance to the study 
of abnormalities phenomena as well as the design process. The more accurate the model, the more reliable 
and accurate the results [12]. Various researchers use optimization approaches to solve many problems by 
estimating parameters [13]-[16]. Many techniques are used to improve the estimation of the transformer's 
parameters by employing various optimization strategies for estimating the parameters from search space, 
constraints, and objective functions [17]-[19]. The estimated coyote optimization algorithm (COA) 
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parameters lead to a high closeness to the experiments that achieve the most effective parameters compared 
to other optimizing algorithms between the parameters estimated and the actual parameters [17]. By 
minimizing certain objective functions, particle swarm optimization, and genetic algorithm are used to track 
nameplate data. The findings show that transformer equivalent circuit parameters can be precisely defined by 
evolutionary computation techniques [18]. The imperialist competitive algorithm (ICA) and gravitational 
search algorithm (GSA) are proposed for estimating transformer parameters from nameplate data. These 
techniques can give better accuracy in estimating the parameters of power transformers [19]. The load test 
data is adequate to estimate the parameters using a genetic algorithm at a specific operating point: voltage, 
load currents, input power, and load impedance. The process is used for the determination of output, power 
system load flow, and design of three winding Transformer protective circuits [20]. The paper describes the 
method for determining the parameters of the non-linear current-flow relation which characterizes the 
saturation. They suit very well with the analytical role for modifying the experimental measurements [21]. To 
optimize the transformer parameter, the proposed algorithms use manta ray optimization and chaotic manta 
ray. In contrast with the approach suggested in the previous article, the proposed objective function does not 
require any loss in the estimation process. The approximate parameters obtained through the proposed 
optimization and objective function are compared with the values obtained using the classical Institute of 
Electrical and Electronics Engineers (IEEE) test protocol and the values derived using the literature 
previously introduced. The measurements of the experimentally defined transformer parameters shall be 
performed to demonstrate the effects of the proposed parameter estimation algorithm and the correlation shall 
be studied [22]. The proposed turbulent flow of water-based optimization (TFWO) algorithm used to 
estimate the transformer parameters of the power transformer from its nameplate data and the accuracy and 
reliability of the proposed algorithm are higher than the other algorithms. 

The following is how this paper is structured; section 2 explains the analysis of objective function 
that solves problem formulation. Section 3 contains the details of the proposed TFWO algorithms. Section 4 
analyses the results of the case study and the conclusion is in section 5. 


2. PROBLEM FORMULATION AND OBJECTIVE FUNCTION 
2.1. Transformer model 

The equivalent circuit referred to the primary side of a single-phase two-winding transformer is 
shown in Figure 1. The conventional method of estimation of transformer parameters requires usually open 
circuit and short circuit test experimental data. In this paper, the challenge is to estimate transformer 
parameters by only using the nameplate data without performing tests. The proposed approach is based on a 
reduction in the square error of the nominal transformer nameplate parameters and the corresponding 
parameters calculated in Figure 1. Use the circuit laws of Kirchhoff by solving several mathematical 
equations [23]. The total impedance of transformer is as follows: 


= ZoX(Z2+Rioad) 
ee a = Zo+Z2+Rioad o) 
where: 
Z, =R, + jX, Ź = Á, + jXp, and Z, = on 
from Figure 1, the primary current is: 
za 
L= 2) 
and the secondary current is: 
NR Zo vı 
i ~ Zo+Z2+Rioad 2 Z @) 
the secondary voltage calculated as: 
Vz = Í X Rioaa (4) 
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Where: 

R, = Primary winding resistance 

R, = Secondary winding resistance referred to primary side 
X, = Primary leakage reactance 

X, = Secondary leakage reactance referred to primary side 
R, = Resistance corresponding to core losses 

Xm = Magnetizing reactance 

I, = Primary current 

Í, = Secondary current referred to primary side 

V, = Terminal voltage on the primary side 

Ý, = Secondary terminal voltage referred to primary side 
I, = No — load current 

Rioaa = load resistance 


Figure 1. Referred to the primary side transformer equivalent circuit [14] 


2.2. The objective function of estimation 

The optimization algorithms are used in extracting the parameters of the different ratings of the 
transformer (4 kVA, 10 KVA, and 15 kVA), these techniques need an objective function for performing the 
optimization process. Minimizing the sum of square error (SSE) is the objective function of the optimizer 
technique. The SSE function includes the summation of the square error for the primary current and 
secondary current and voltage referring to the primary. The mathematical formula to compute SSE is as 
follow: 


SSE = Ch act = Best)” t (D act = Dest)? + (Vzact B Vzest )” (5) 


where the first term in the SSE is the difference between the primary actual and estimated current, the second 
term is the difference between the actual and estimated secondary current, and the third term is the difference 
between the secondary actual and estimated voltage. The optimization algorithm needs boundaries limit; 
these boundaries are explained in Tables 1-3 [22], [23]. 


Table 1. The extracted parameters boundaries for Table 2. The extracted parameters boundaries for 
4 KVA 10 kVA 
Parameters Lower bound _ Upper bound Parameters Lower bound _ Upper bound 
Ry 0.2 0.6 Ry 0.3 1 
xX, 0.1 0.3 xX, 0.2 1 
R,' 0.2 0.45 R,' 0.1 1 
Xo! 1.5 2.5 Xx, 0.2 1 
Re 1400 1500 Re 500 1000 
Xm 700 750 Xm 200 500 


Table 3. The extracted parameters boundaries for 15 kVA 


Parameters Lower bound Upper bound 
Rı 1:2 2.5 
Xı 1.3 3.5 
R,' 1.3 2.2 
X,' 1.8 2.3 
Re 99,000 12,0000 
Xm 9,000 9,200 
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3. TFWO ANALYSIS 

A circular shape is formed from the turbulent flow of water; this circular form is called a whirlpool. 
gravity force is affected on this whirlpool so that water forms a spiral path [24]. The random nature behavior 
in oceans, rivers, and seas is the feature of the whirlpool phenomenon. The suckhole that sucks the objects in 
the center of the whirlpool. At first, the initial population (X°, comprising N, members) of the algorithm is 
divided equally between Nwpn groups or whirlpool sets, and then the strongest member of each whirlpool set 
(The better objective value member f (X) is seen as an object pull whirlpool (X, including Np - Nwn objects). 

Every whirlpool (Wh) behaves as a sucking well, and tends to unify the locations of objects inside 
its set (X) with its central position by applying a centripetal force on them, and push them into its well. Thus, 
the jth whirlpool and the local position on W hj combine the ith object position (X;) with itself ( X; = Wh, ). 
However, other whirlpools produce some deviations (AX;) because of the distance among them (Wh — Wh,) 
and its objective values (f(X)) as well. Accordingly, the new position of the ith object becomes X/’” = 
Wh, — AX;. and the objects (X) travel through their whirlpool's core and toward it with their unique angle 
(ô). Hence, this angle at each iteration is changing according to (6): 


oy’ = 6; + rand, * rand, * T (6) 
to model and calculate the farthest and nearest whirlpools (AX; ), in (7) depicts the whirlpools with the least 


weighed distance from all objects, and then AX; is calculated using both in (8). In (9) is used to update the 
position of the particle. 


A= f(Wh,) * |Wh, — sum(X;)|°° (7) 
AX; = cos(6;"*”) * rand(1. D) * (Wh; — Xi) — sin(6?*”) * rand (1. D) * (Wh, — X;)) 

* (1 + |cos(67°") — sind?) |) (8) 
KOF = Wh, _ AX; (9) 


where Why and Why manifest the whirlpools with minimum and maximum of A+, respectively, while 
6; characterizes the ith object’s angle. Centrifugal force (FE;) sometimes overcomes the centripetal of the 
whirlpool, and randomly transfers the object to a new location. The centrifugal force is modeled as illustrated 
in (10), which randomly occurs in one dimension of the decision variables. The centrifugal force is 
determined based on the angle between the whirlpool and the object, as shown in (10), and if this force is 
greater than a random value in the range [0, 1], the centrifugal action is performed for randomly selected 
dimension, as shown in (11). This phenomenon can be expressed mathematically as follows: 


FE; = ((cos(67™))* * (sin(67™))?)? (10) 

Xip = Xp" + rand » (Xp — xiv) (11) 
the whirlpools collide and displace one another. This phenomenon can be modeled as the same as the impacts 
of whirlpools on the objects, where every whirlpool tends to pull other whirlpools and imply the centripetal 


force on them. As shown in (12), the closest whirlpool can be mathematically represented using the minimum 
sum and its objective function. Then, the whirlpool’s position can be updated according to (13) and (14): 


A= f (Wh,) * |Wh, — sum(Wh,)| (12) 
AWh, = rand(1.D) + |cos(7™) + sin(67™)| * (Wh; = Wh;) (13) 
Whee = Why a AWh,; (14) 


where 6; represents the jth whirlpool holes angle value. Finally, the strongest member of the new members 
of the whirlpool's set is chosen as the new whirlpool for the next iteration when the value of the objective 
function is less than the preceding whirlpool. The TFWO flowchart is depicted in Figure 2. 
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Figure 2. Flowchart of the proposed TFWO algorithm 


4. ANALYSIS OF RESULTS 

This section presents an analysis of the parameters extracted using the proposed TFWO algorithm 
for various transformers rating. The objective function is applied in this paper on the main proposed 
algorithm TFWO and it's applied also on the Cuckoo search algorithm (CSA) [25]. The performance of 
TFWO is validated by comparison of its results with this CSA and other algorithms from the previous work 
such as genetic algorithm (GA), practical swarm optimization (PSO) [18], imperialist competitive algorithm 
(ICA), and gravitational search algorithm (GSA) [19]. 
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4.1. Results for transformer types 

In this subsection, the details of extracted parameters and the best objective function by TFWO is 
performed. The proposed TFWO algorithm is compared with the cuckoo search algorithm (CSA) for 30 
independent runs. The proposed algorithm test is applied to transformers rated at 4 kVA, 10 kVA, and 
15 kVA, and the results are compared to existing approaches. 


4.1.1. The transformer of rating 4kVA 

Parameters were extracted from a single-phase transformer based on the analysis in section two and 
three. Table 4 explains the value of the parameters for the 4 kVA transformer. Table 5 explains the current 14 
and I’, and the voltage V>» for this transformer based on the parameters extracted from each algorithm at the 
best objective function. The value of the best SSE for each algorithm is illustrated in Table 6. The statistical 
results of 30 run for TFWO and CSA algorithms are presented in Table 7. From the results in Tables 4-6, it 
can be observed that the proposed algorithm TFWO is the best algorithm in extract transformer parameters 
and the CSA algorithm is the second algorithm then ICA, GSA, PSO, and GA respectively. From Table 7 the 
accuracy and reliability of the TFWO algorithm are higher than the CSA algorithm. 


Table 4. The parameters extracted for a4 kVA transformer at the best SSE 


Algorithm actual TFWO CSA GA PSO ICA GSA 
R1 0.4 0.387421895 0.35306923 0.598 0.587 0.430 0.425 
X1 0.2 0.1 0.298748267 0.226 0.2554 0.202 0.203 
R2' 0.4 0.45 0.45 0.336 0.209 0.394 0.415 
X2' 2 25: 2.5, 1.957 1.602 2.5 2.399 
Re 1500 1400 1400 1410 1476 1200 1426 
Xm 700 700 700 707 738 700 745.3 


Table 5. The value of the current and voltage for the 4 kVA transformer of each algorithm 


Transformer data Actual TFWO CSA GA PSO ICA GSA 
Il 14.0813 13.91535929 13.91535929 14.1035 14.0818 13.82729 13.8077 
12' 13.6893 13.68986622 13.68986622 13.6654 13.6972 13.6196 13.5916 
V2' 235.8759 235.8856563 235.8856563 234.2131 237.2283 235.8917 235.407 


Table 6. The value of the best objective function for 4 kVA transformer of each algorithm 


Algorithm The best SSE 
TFWO 0.027631826 
CSA 0.027631826 
GA 2.76596789 
PSO 1.82904842 
ICA 0.06962881 
GSA 0.30426946 


Table 7. Statistical analysis of fitness function for 4 kVA transformer 


Objective function (SSE) TFWO CSA 
Minimum 0.027631826 0.027631826 
Mean 0.027631826 0.027631826 
Maximum 0.027631826 0.02763 1827 
Standard deviation 2.31E-15 2.19E-10 


4.1.2. The transformer of rating 10 kVA 

One phase transformer of rating 10 kVA, 50 Hz, 500/125 V is used to extract their parameters based 
on the analysis in section two and three. Table 8 explains the value of the parameters for the 10 kVA 
transformer. Table 9 explains the current J, and l’, and the voltage V» for this transformer based on the 
parameters extracted from each algorithm at the best objective function. The value of the best SSE for each 
algorithm is illustrated in Table 10. The statistical results of 30 run for TFWO and CSA algorithms are 
presented in Table 11. From the results in Tables 8-10, it can be observed that the proposed algorithm TFWO 
is the best algorithm in extract transformer parameters and the CSA algorithm is the second algorithm then 
PSO, GA, ICA, and GSA respectively. From Table 11 the accuracy and reliability of the TFWO algorithm 
are higher than the CSA algorithm. 
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Table 8. The parameters extracted for 10 kVA transformer at the best SSE 


Algorithm Actual TFWO CSA GA PSO ICA GSA 
R1 0.9 0.3 0.300000525 1.025 0.811 0.8 0.8001 
X1 0.94 0.315511728 0.31550416 0.8 0.8608 0.8 0.8119 
R2' 1.6 0.1 0.1 1.507 1.678 1.5 1.5004 
X2' 0.44 0.2 0.2 0.493 0.7540 0.4259 0.4236 
Re 700 1000 1000 651.5 713 692.48 695.54 
Xm 250 500 500 204.4 314.2 255 251.35 


Table 9. The value of the current and voltage for 10 kVA transformer of each algorithm 


Transformer data Actual TFWO CSA GA PSO ICA 
Il 19.688 19.82456533 19.82456533 19.6072 19.7012 19.0427 19.0417 
12' 19.299 19.29878424 19.29878424 19.1979 19.3189 18.2217 18.2207 
V2' 491.729 491.7235026 491.7235026 487.847 487.9916 455.5431 455.5186 


Table 10. The value of the best objective function for the 10 kVA transformer of each algorithm 


Algorithm The best SSE 
TFWO 0.018680358 
CSA 0.018680358 
GA 15.08667385 
PSO 13.96872901 
ICA 1310.996346 
GSA 1312.773503 


Table 11. Statistical analysis of objective function for 10 kVA transformer 


Objective function (SSE) TFWO CSA 
Minimum 0.018680358 0.018680358 
Mean 0.018680358 0.018680358 
Maximum 0.018680358 0.018680358 
Standard deviation 1.84E-16 4.99E-15 


4.1.3. The transformer of rating 15 kVA 


645 


One phase transformer of rating 15 kVA, 50 Hz, 2400/240 V is used to extract their parameters 
based on the analysis in section two and three. Table 12 explains the value of the parameters for the 15 kVA 
transformer. Table 13 explains the current J} and l’, and the voltage V» for this transformer based on the 
parameters extracted from each algorithm at the best objective function. The value of the best SSE for each 
algorithm is illustrated in Table 14. The statistical results of 30 run for TFWO and CSA algorithms are 
presented in Table 15. From the results in Tables 12-14, it can be observed that the proposed algorithm 
TFWO is the best algorithm in extract transformer parameters and the CSA algorithm is the second algorithm 
then PSO, GA, GSA, and ICA respectively. From Table 15 the accuracy and reliability of the TFWO 


algorithm are higher than the CSA algorithm. 


Table 12. The parameters extracted for a 15 kVA transformer at the best SSE 


Algorithm Actual TFWO CSA GA PSO ICA GSA 
Rl 2.45 1.2 1.216390598 2.76 2.25 2 2 
X1 3.14 2.130401219 1.811673414 3.414 4.082 3 3.11 
R2' 2 1.3 1.3 1.68 2.2 1.8 1.81 
X2' 2.2294 1.8 1.8 1.846 1.8526 2 2.26 
Re 10,500 12,000 12,000 97,001 99517 120000 104281 
Xm 9,106 9,200 9,200 8,951 9,009 9,200 9094.87 


Table 13. The value of the current and voltage for the 15 kVA transformer of each algorithm 


Transformer data Actual TFWO CSA GA PSO ICA GSA 
Il 6.2 6.226571224  6.226571224 6.2017 6.2004 6.1653 6.1693 
2' 6.2 6.199999819  6.199999819 6.2001 6.2008 6.1387 6.1393 
Vv2' 2383.8 2383.799931 2383.799931 2381.8 2384.7 2375.662 2375.917 
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Table 14. The value of the best objective function for the 15 kVA transformer of each algorithm 


Algorithm The best SSE 
TFWO 0.000706035 
CSA 0.000706035 
GA 4.0000029 
PSO 0.8100008 
ICA 66.23200578 
GSA 62.14631598 


Table 15. Statistical analysis of objective function for 15 kVA transformer 


Objective Function (SSE) TFWO CSA 
Minimum 0.000706035 0.000706035 
Mean 0.000706035 0.000706035 
Maximum 0.000706035 0.000706035 
Standard deviation 1.00E-17 1.31E-12 


5. CONCLUSION 

This paper proposes a turbulent flow of water optimization for transformer parameter extraction 
using nameplate data, which is applicable to single-phase and three-phase transformers. The parameters 
extracted from the proposed TFWO algorithm are executed successfully for three separate transformer 
ratings. The results of comparisons between TFWO, ICA, GSA, PSO, and GA estimate errors in the 
transformer parameters show that the TFWO results are more accurate and reliable than the ICA, GSA, PSO, 
and GA estimate errors. The transformer nameplate and the search space direct the TFWO method's objective 
function, ensuring that the proposed method is accurate, reliable, and complete. As compared to other 
methods, the TFWO method has the lowest objective function standard deviation of summation of square 
errors for primary current and secondary current and voltage corresponding to the primary. If experimental 
testing is not possible, the proposed TFWO technique can execute for the transformer three times instead. 
The proposed TFWO technique can be adjusted for more than six parameters optimization. Then for the 
three-phase transformer model includes distributed resistance and inductance, mutual inductance and three 
capacitances including self and ground capacitance and inter-winding may be optimized in future work. The 
results of the proposed method are compared with actual measurements obtained from previous publications 
to prove the validity of the proposal. 
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